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We improve the estimate of the axion CDM energy density by considering the new values of 
current quark masses, the QCD phase transition effect and a possible anharmonic effect. 

PACS numbers: 14.80.Mz, 12.38.Aw, 95.35.+d 



OO 
O 
O 

(N 
C 

On 



Or 
i ' 



> 

o 

\6 
o 

OO 

o 



X 



I. INTRODUCTION 

The standard model is the greatest triumph of the 20th 
century particle physics. However, it suffers several nat- 
uralness problems. The so-called strong CP problem is 
one of them. Up to now, the most promising strong CP 
solution is axion [l[ . Although its coupling is very small, 
it could have affected the evolution of the Universe. 

On the other hand, the observations like the galaxy 
rotation curve and the cosmic microwave background ra- 
diation tell us that a significant portion of the energy den- 
sity of the universe cannot be explained by the ordinary 
baryonic matter. The recent analysis of WMAP data[2| 
combining the observation of the large scale structure 
prefers the standard vacuum energy and cold dark mat- 
ter (ACDM) model, with the cold dark matter (CDM) 
density fl DM h 2 ~ O.f f 43 ± 0.0034. 

During the inflation, a sufficiently light scalar field be- 
comes spatially homogeneous and its value is selected 
stochastically Q. After the inflation, the Hubble param- 
eter decreases as the universe expands. If the Hubble 
parameter becomes smaller than some scale, usually the 
scalar mass scale, the scalar field starts to roll down. In 
the hydrodynamic description, this coherent oscillation 
works as cold matter in the evolution of the universe. If 
the interaction is small enough to keep coherence until 
the matter-dominant era, it can be a part of the present 
matter energy density of the universe. The axion is that 
kind of particle and can be coherent until even now. 

Just after the invisible axion was introduced Q , people 
recognized its cosmological implications that the classical 
axion field oscillation behaves like CDM [f| . On the axion 
CDM energy density, Turner presented a more accurate 
number including anharmonic effect and numerical in- 
stanton calculation in late 1980s DeGrand et. al. 
concerned about the possible effects of the chiral phase 
transition Q- 

After Turner's pioneering work, many input parame- 
ters in axion cosmology, especially the QCD parameters, 
are refined. We present the improved number for the ax- 
ion dark matter relic density with these new parameters. 



We also discuss the possible effects of the chiral phase 
transition on the total entropy and the axion number 
density. This occurs only when supercooling arises. 



II. COMPUTATION OF AXION RELIC 
DENSITY 

Axion is a pseudo Nambu-Goldstone boson, acquiring 
mass from the instanton effect. The temperature depen- 
dence of the axion mass is important in calculating axion 
energy density for which we will use the result of Ref. § 
obtained in the dilute gas approximation. For the pur- 
pose of sketching the computation, we just parameterize 
the axion potential as 



V{6) = -C(T) cos(0), 
where 8 = a/ F a and 

C{T) = a inst GeV\T Gc y)~ n 



(f) 



(2) 
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where T GoV = T/GcV. 

The inflation freezes all fields whose masses are smaller 
than the Hubble parameter at that time. It also flattens 
the fluctuation of the field. Consequently, just after the 
inflation, the sufficiently light fields remain homogeneous 
with a certain value that is determined stochastically. 
The axion and the moduli are typical representatives of 
them. 

At the end of inflation, the universe is reheated by some 
processes. Since axion is very weakly interacting, the re- 
heated radiation of the early universe cannot destroy its 
coherence. As the universe cools down, the Hubble fric- 
tion can no longer be the dominant term and the axion 
field starts rolling down the potential hill and continues 
to oscillate. Since the Hubble parameter becomes smaller 
and smaller than the axion mass, or the oscillating period 
of the axion is smaller than the Hubble time, its energy 
and pressure can be described by the time averages dur- 
ing each oscillation. In the case of the harmonic poten- 
tial, the effective pressure is zero and the axion field acts 
like cold matter in the evolving universe. Moreover, if the 
axion can keep the coherence until the matter domina- 
tion epoch, it can contribute to the present CDM energy 
density. 

Since the axion mass depends on the temperature or 
the time, the comoving energy density conservation does 
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not work in estimating the axion relic density. How- 
ever, the adiabatic theorem says that as long as the 
adiabatic conditions H,m a /m a <C m a hold, we can 
find the adiabatic invariant I. For the harmonic poten- 
tial, the invariant is the comoving axion number den- 
sity / = R 3 p a /m a (T). Since the coherent state is not 
the number eigenstate, the interpretation as the number 
density only works for the harmonic potential, i.e. the 
free field case. We will discuss the appropriate invariant 
quantity in the case of the interacting axion field. 



III. TEMPERATURE DEPENDENCE OF 
AXION MASS 



The QCD vacuum has the quark-gluon phase and the 
chiral symmetry breaking hadronic phase, separated by 
the critical temperature T c of the chiral symmetry break- 
ing. T c has been calculated to be 148±|?(172±|2) MeV 
for three (two) light quarks @- This region is in the 
boundary of weak and strong coupling regimes and it is 
difficult to estimate the axion mass very accurately in 
this region. Early attempts to estimate the temperature 
dependent axion mass has been given in 0, [ll[ . In these 
earlier studies, this phase transition was not considered 
and here we point out that the consideration of the phase 
transition docs not lead to a substantial change in the es- 
timation of the axion energy density. 
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FIG. 1: The phase transition near the critical temperature 
T c ~ 150 MeV. The solid line describes a smooth path first 
order transition and the dashed line describes a supercooling 
case. 

In the hadronic phase below the critical temperature, 
the correlation length of q — q pair is large and we assume 
a temperature independent height for the axion poten- 
tial as schematically depicted in Fig. [T] in the sudden 
phase change approximation. In this phase, the squared 
axion mass has the well-known one power of the light 
quark mass, and the coefficient of the axion potential is 
fl<m%Zl{\ + Z) 2 [ll which is the second line in Fig. [2]m 
the expansion of the determinental interaction [131 ]. On 
the other hand, in the quark-gluon phase mesons such 
as 7T° and rf are not dynamical, and the first term of 
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FIG. 2: The determinental interaction of light quarks. 



Fig. [2] is the coefficient of the axion potential, which is 
—K~ 5 uuddss cos 9 for three light quarks where 9 is the 
QCD vacuum angle. Considering the 3x3 mass ma- 
trix for a,7r ,?/, the insertions of qq VEVs (the lightly 
shaded diagram of Fig. [2]) is the dominant contribution 
to the n' mass [lH from which we estimate the strength 
K~ 5 w m 2 ,f 2 , /4v 9 ~ 10- 13 McV" 5 for f v , ~ / w in the 
hadronic phase and v f=a 260 MeV. The same number A'~ 5 
also appears in the quark-gluon phase. This value will be 
needed if the adiabatic approximation is violated, but we 
will show below that the adiabatic approximation is not 
violated in the bulk of the remaining parameter space 
and the number K will not be used in the end. 

Closing the quark lines to insert the current 
quark masses in the quark gluon phase, we obtain 
—K~ 5 (m u mdm s /(ft) cos 9 where p is the effective instan- 
ton size in the instanton size integration at high temper- 
ature. Instead of this naive averaging, we use the original 
expression given in Q 

Tl(p, 0) e {-^ 2 (™+N f )-12A(\)[l+±(N-N f )}} (3) 

where A = npT,A(X) ~ -^ln(l + A 2 /3) + a(l + 
7A -2/3)-8 with a = 0.01289764 and 7 = 0.15858. The 
zero temperature density is 



n(p,0) =m u m d m a C N ^pfj s (±£) 2N e -^^ 2 ^ 



(4) 
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with £ = 1.33876 and C N = 0.097163 for N = 3. We 
use A = 2.5 GeV and N = 3 where a c (A = 2.5 GeV) ~ 
0.28 ±0.01 



a c (/i) 



4tt 



■ - /3oln( M 2 /A^ CD ) 



2/3i ln[ln(Ai 2 /A^ CD )] 
#T '"(A'VA^cd) 



In 



~8£f 



ln(/z 2 /A 2 



QCD; 



where /3 = 11 - |iV/,/3i = 51 



f JV/, and/3 2 



(5) 



2857- 



^p-Nf + ^§fNj. The leading factors are m u m d m s T and 

(A/T)^ 11N ~ 2Nl \ So, we parametrize the instanton size 
integration of © at T = TqoV GeV near 1 GeV (from 
700 McV to 1.3 GeV) as 

V{6) = -C(T) cos(0), (6) 

where 9 = a/F a and we parametrize C(T) as 

C(T) - « inst GcV 4 (T GcV )-". (7) 

For A QCD = 380 (440, 320) MeV 0, we obtain 
a inst = 3.964 x 10- 12 (1.274 x 10" 11 , 9.967 x 10~ 13 ), and 
n = 6.878 (6.789, 6.967). It starts rolling down the axion 
potential hill when the axion mass term becomes dom- 
inant over the Hubble friction term in the equation of 
motion, i.e. 3H m m(T). The Hubble parameter is 
lMgl /2 {T)T 2 /M P . At T w lGcV, u, d, s, e, fi, 2,uD, 8G, 7 
remain relativistic. So g*(T w lGeV) = (8 + 1) x 2 + 
|{(3 x3 + l + l)x2 + 3}x2 = 61.75. Equating 3H{T) 
and m(T) = y/C(T)/F*, we obtain 



Ti,GeV — < 



0.808 x 



0.916 x 



1.020 x 



F a , G c 



Fg, Go 



-0.184 



-0.185 



(8) 



for Aqcd 
tively. 



320 McV, 380 McV, and 440 McV, respec- 



IV. ANHARMONIC EFFECT ON AXION 
NUMBER 

In early studies [B[ , the anharmonic higher order terms 
in the axion potential were neglected and only the har- 
monic term (the mass term) was considered. In such 
an approximation, the axion total number in the comov- 
ing volume is conserved whenever the adiabaticity con- 
ditions, H,m a /m a <C m a are satisfied. Turner @, |Tl| 
considered the effect due to the anharmonic terms on 
the axion CDM energy density. Later, Lyth presented 
an extensive study on this issue [l5[ . In the early epoch 
after the classical axion field starts to roll down, the axion 



number is not conserved because the anharmonic terms 
are not negligible when 9i > 1. Physically, the introduc- 
tion of the anharmonic terms mean the axion number 
changing interactions so that the axion number conser- 
vation is not guaranteed. 



A. Adiabatic invariant 

Here, we present the method to treat this anharmonic 
effect. The axion field satisfies 



3H9 + V = 



(9) 



where 9 = a/F a , V(9) = V(9)/F 2 , and the prime denotes 
the derivative with respect to 9. This equation of motion 
can be derived from the time-varying Lagrangian, 



L = R 3 (^F 2 9 2 -V 



(10) 



where V = m 2 F 2 (l — cos#). We treat here R and m a as 
time- varying parameters. Under the adiabatic condition, 
H,m a /m a <g; m , we have the adiabatically invariant 
quantity [l6|, 



(ii) 



where q = 9 and p = R 3 9 which is a conjugate momen- 
tum of q. Using this quantity, we can derive an invariant 
quantity of the axion evolution, 



where 



A(0) 



R 3 m a 9 2 f 1 (9) = constant 
2V2 



(12) 



Vcos 9' - cos 9 d6'. (13) 



The numerical result is shown in Fig. [3] The same form 
was also obtained in flal. 




2.0 2.5 3.0 
FIG. 3: The correction factor from anharmonic effect. 
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Thus, wc obtain 



p a (T 7 ) ~ m a (T 7 )m a (Ti) 



( R 3 (Ti) 
U 3 (T 7 ) 



%h(0i), (14) 



where T 7 = 2.73K. We approximate here /(0(T 7 )) ~ 1 
since axion amplitude decreases as the universe expands 
and present axion amplitude is small enough to be har- 




Initial anharmonic correction 



There exists the adjustment in the initial stage of the 
rolling down. This initial anharmonic correction will turn 
out to be larger than the general correction presented in 
the previous subsection. Wc assumed two adiabaticity 
conditions. m a /m a <C m a and H <C m a . However, since 
the expansion rate during the first few oscillations are 
not so small, 3 if ~ m a , we should be careful in treating 
the anharmonic correction in the initial stage. Actually, 
even for the case of the harmonic limit, 9\ <C 1, the 
initial adjustment is needed. For a large initial misalign- 
ment angle, 9\ ~ 1, the duration time during which the 
expansion rate is not so small is longer than for the case 
of a small 6q . The effect of the initial adjustment is larger 
for the larger 9\ . 

To estimate the initial adjustment correction, we solve 
the equation of motion numerically and obtain another 
temperature T 2 and the corresponding misalignment an- 
gle 62 to use them in the subsequent calculation rather 
than using the original T\ and 9\ . We have noted that a 
consideration of the first half oscillation is enough for this 
purpose toward obtaining H <g; m a as depicted in Fig. 
2J The common tangent (the dashed line) to the dashed 
curves is the square root of the invariant y/l. In fact, this 
invariant shows an overshoot (a factor of 1.8 which will 
be shown later) from the initial misalignment angle 9\, 
the initial value of the solid curve, which occurs because 
in this initial stage the expansion rate is non-negligible 
compared to m a . Hence we obtain a larger axion CDM 
energy compared to the previous estimates. 

For the presentation of the result, wc parametrize 
h/ti = / 2 ~ 2 (6'i,n) and 62/81 = fz{9\,n) where n is 
the exponent in the temperature dependent axion mass 
discussed earlier. The correction factors / 2 (#i,n) and 
fz{9\,ri) are shown in Fig. \E\ 

If the significant entropy production is absent dur- 
ing the cooling period from T\ down to T 2 , T2/T1 = 
(WO = f2{fii,n). To split the pure anhar- 
monic effect from the rest, we parameterize these as 
T2/T1 = c 2 / 2 (0i, n) and 6 2 /0i = c 3 f 3 (9i, n) where c 2 = 
lim 01 ^ o h — 0-691 and c 3 = lim^o f 3 ~ 0.414. 

Then, the appropriate 6% dependence on the relic den- 




m a (T l )t 



FIG. 4: The solid curves represent the 6 evolution for 8\ =0.1 
(in the upper figure) and 81 = 3 (in the lower figure). The 
common tangents to the dashed curves are the square root of 
the invariant \fl. After a half period, they almost approach 
to the condition H <g m a - 



sity is given by 

fia « ^3 /i(0 2 ) 

9j4f§(Sh, n)/ 1 (0 1 c 3 /3(gi,n)) 
4 +n/2 f 3 2 +n/2 (0 u n) 
= 1.846 x 9l x F(9 l ,n) 



(15) 



where the function F (6q,n) = 

fl{9i,n)h(0 1 c 3 f 3 (0 1 ,n))/fl +n/2 {9 1 ,n) which is shown 
in Fig. [6] 



V. QCD PHASE TRANSITION AND AXION 
NUMBER CHANGE 



We know that the adiabatic invariant I is conserved in 
the comoving volume. We calculated this at a later time 
t 2 > ti where rn a /m a <C m a are satisfied. Thus, if 
we know the initial misalignment angle 9\ = 9{T\), we 
can determine the current axion energy density. Dur- 
ing the universe evolution from T = T\ to T = 2.73 K, 
however, there is the QCD chiral phase transition epoch 
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FIG. 5: The initial correction factors /2 and 
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FIG. 6: The combined correction factor F(8i,n) 



when quarks and gluons in the quark-gluon phase com- 
bine each other and change to hadrons in the hadronic 
phase. Hence, we must check the effect of this phase 
transition on the total axion number. Since we do not 
know the exact phenomena during this phase transition, 
we must interpolate two regions of the quark-gluon phase 
and the hadronic phase as shown in Fig. [TJ 



A. Possible supercooling effect in phase transition 

In 0, 03 > DeGrand et al. studied supercooling ef- 
fect during the QCD phase transition for various cases. 
We briefly explain their argument. Assuming that the 
universe undergoes first-order phase transition at T ~ 
Aqcd, we shall follow bag model of phase transition. In 
the bag model, the equation of state is given by Fig. [7] 



p q (T)/B=[3r/(r-l)]T 4 + l, 
p q (T)/B=[r/(r-l)]f\ 
T cSq (T)/B = [4r/(r - l)]f 3 , 



in the quark-gluon phase and 

Ph (T)/B = 3p h (T)/B 
T c s h (T)/B = [4/(r 



[3/(r 
Dlf 3 , 



l)}f 4 



(16) 
(17) 
(18) 

(19) 
(20) 



in the hadronic phase, where T = T/T c , r = g q / gh — 3, 
B = (gh — g q )(7r 2 /90)T ( f and gq(<7h) is the relativistic de- 
grees of freedom of quark-gluon(hadron) phase. We shall 




Tj 
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(b) 



FIG. 7: The first order phase transition. 
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interpolate p and p linearly between 7] and Tf like Fig. 
[7J Subscript i(f) denotes the initial(final) temperature 
during the phase transition. After some simple algebra, 
we obtain the following relations [13] 



Ri/R, = [(pi - B)/ Pf ] 



(Pf-Pi)/4(p f -Pi + B) 



x(t{-t i ) = -VB(p { -p i ) 



dx 



Bp t - (pi - Pi + B)x 2 ' 



SSI Si = r-^ 4 [( Pi - B)/p f } 3B / i ^- B -P'1 

where Pi = p q (Tf), p f = p h (T f ), and X = HitGnB) 1 / 2 » 
0.22 GeV 2 /TOpi- Since we do not have enough informa- 
tion on Ti and Tf, we treat these as free parameters. 
Then, we obtain the phase diagram in the T;-Tf plane 
as shown in Fig. [5] 




FIG. 8: The physically allowed region in the TJ-Tf space. 
Dashed lines (with numbers shown below) stand for contour 
of xAt, and the red wedge for an exaggerated possible sudden 
change is x^t < 10 -3 . The violet shaded regions show the 
entropy production As/s; during the phase transition. 



B. Axion number change 

Now we estimate axion number change in the allowed 
region. To do this, we shall examine if the adiabaticity 
condition is still valid or not in this region. Since H is 
decreasing during the expansion of the universe, m a > H 
condition is satisfied after the axion field starts to roll 
down. We must check whether the adiabaticity condition 
ni a /m a < m a is satisfied during the phase transition. 



1. Smooth transtion 

For a smooth transition, T\ = Tf = T c , we have p Y 
[3r/(r - 1) + l]B, p f = 3B/(r - 1). Thus, we obtain 



RfJ R[ 



,1/3 



(21) 



X At = -(r - l) 1/2 [arctan(4r - 1) 1/2 - arctanV3] w 0.22, 

(22) 

AS = 0, (23) 

where At = tf — t[. Therefore, we obtain At ~ 
10 10 eV _1 S> m^ 1 . Thus, the adiabaticity condition is 
satisfied and the axion number in the comoving volume 
is not changed. 



2. Abrupt transition 

Along the red wedge in Fig. [HJ pf = pf. Then, 

Rf = Ri, (24) 
At = 0, (25) 
AS /Si = Tf/rT? - 1. (26) 

In this case, At — and hence we cannot sustain the 
adiabaticity condition. Thus, the axion number must be 
changed. Depending on the phase of the axion field, the 
axion number may be either lowered or raised by factor 
0.1 or 10 0- 



Parameters in colored region are not allowed because 
of some physical constraints. The dark blue region is 
excluded from the condition that Tf is equal or larger 
than T; . The green region is excluded from the condition 
that entropy never decreases. The yellow region(beyond 
the red line) is excluded from the constraint that energy 
conservation is not violated. The light blue region is 
excluded from the constraint on the kinematical limit 
of the bubble nucleation rate[l7j. The thin red wedge 
for X At < 10~ 3 is the exaggerated region (an abrupt 
change case is for X At < 10~ 4 which is not noticeable in 
the figure) where the adiabaticity condition is violated, 
which is explained below. 



3. Bulk region 

In the bulk region in Fig. [51 we show a reliable re- 
gion satisfying the adiabaticity condition. In Fig. [8j we 
can see that very tiny part (red wedge) of the bulk re- 
gion breaks the adiabaticity condition. Thus, the axion 
number is not changed in the allowed region except the 
red wedge. In Fig. [8j the entropy increase is shown as 
violet color shades. For a small(< 50%) entropy produc- 
tion during the transition in this region, the axion energy 
density changes because the entropy increase ratio 7 is 
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changed. Entropy change in the bulk region is given by 



As 

Si 



-1/4 



rpA 



- 1. 



(27) 



Thus, we conclude that the chiral phase transition does 
not have any serious effect on the initial total axion num- 
ber in the comoving volume. 



VI. RESULT 



The axion energy density, using the new numbers on 
the current quark masses [181 ] and including anharmonic 
effect and QCD phase transition effect, is estimated as 



J 



p (T 7 = 2.73K) = TO a (T 7 K(T 7 )/i(# 2 ) = 



l + Z 



Z , g.sjTjT* 3 x 1.66 8 2 2 h{0 2 ) F t 



7 



Alp VTi 



To 



-3-n/2 



1 



0.785 x 1CT 11 



1.449 x 10~ n 



F a , Ge y\ 1 el fl(e 1 ,n)h(e 1 f 3 (e 1 ,n)) 



10 1 



Ti,GcV 7 



3+n/2 



(On n) 



(28) 



F, 



a,GcV 



1 



10 12 J Ti GeV 7 



lF(e^n)(eV) 4 



where we used Z = m u /md = 0.5, m w o =135.5 MeV, 
=93 MeV and g* s (present)=3.91, is the anhar- 

monic correction given in Eq. (|13p and 7 is the entropy 
increase ratio. Inserting Eqs.© to l]28p. we obtain for 
01 < 1 



(eV)^ 



1.794 x 10- n ^F(6>i) 



1.420 x 10- u ^F(6»i) 



1.580 x lO-n^F^)) 



Fa.GoV 



1.184 




FJy (GeV) 
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FIG. 9: The bound from overclosure of the universe. The 
yellow band shows the error bars of A and two red dashed 
lines are the limits of the allowed current quark masses. The 
anharmonic effect is taken into account, including the initial 
correction factor of Eq. (|15[) . Here, the entropy production 

(n+4)/(n + 6) ^ 



ratio 7 is absorbed into the bracket of F a 
7 - 84 . 



7 = T 



for Aqcd = 320 MeV, 380 MeV, and 440 MeV, re- 
spectively. For the critical density p c = 3.978 x 
10~ n (/i/0.701) 2 (eV) 4 , the axion energy fraction is given 
by 



0.450 



2 / \ 1.182 

0.701 \ I fa.GoV 



0.3971 e -^± 



0.356 



I I To 11 " 
2 , 

0.701 \ / -Fa.GoV 

"7T - ) I To 15 " 

2 / 

0.701 \ I fa.GoV 



(29) 



in 



12 



for Aqcd = 320 MeV, 380 MeV, and 440 MeV, respec- 
tively. From this result, we can determine the cosmo- 
logical bound on the axion decay constant F a which is 
presented in Fig. [9] In Fig. [9l the yellow band shows the 
error bars of A and two red lines are the allowed limits 
of the current quark masses, Z. From the WMAP five 
year data 0, n DM h 2 ~ 0.1143 ± 0.0034, we obtain the 
bound F a < 10 12 GeV for 6 1 ~ 0(1). For 0i/./y < 0.18, 
F a > 10 13 GeV is allowed. The increase of a factor of 
few on the upper bound of F a mainly comes from the 
currently used smaller current quark masses compared 
to the previous ones, and the overshoot of the adiabatic 
invariant \fl from the initial misalignment angle B\ . 
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VII. CONCLUSION 

In conclusion, we presented a new expression for the 
cosmic axion energy density using the new values for the 
current quark masses, implementing the QCD chiral sym- 
metry breaking and the adiabatic invariant /. We in- 
cluded the effect of anharmonic terms also. Using the 
adiabatic invariant quantity /, we presented the correc- 
tion function F to the anharmonic effect and a factor 
1.85 which arose from the overshoot from Q\. Using the 
bag model for the phase transition in the early universe, 
the axion number is shown to be conserved in most cir- 
cumstances. If the axion number change is considered 



during the phase transition, it is extremely fine-tuned in 
the T\ — Tf space, in which case the only change is from 
a small entropy production during the phase transition. 
In this way, we obtain a factor of few larger bound on 
F a compared to the previous ones even for the smaller 
requirement on f2 a . Of course, for a small Q\ a strip with 
large F a s is allowed, which has been used in the anthropic 
arguments [l9| . 
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